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Cross-over in Scaling Laws: 
A Simple Example from Micromagnetics 

Felix Otto* 



Abstract 

Scaling laws for characteristic length scales (in time or in the model param- 
eters) are both experimentally robust and accessible for rigorous analysis. In 
multiscale situations cross-overs between different scaling laws are observed. 
We give a simple example from micromagnetics. In soft ferromagnetic films, 
the geometric character of a wall separating two magnetic domains depends 
on the film thickness. We identify this transition from a Neel wall to an Asym- 
metric Bloch wall by rigorously establishing a cross-over in the specific wall 
energy. 



1. Introduction 

Many continuum systems in materials science display pattern formation. These 
patterns are characterized by one or several length scales. The scaling of these char- 
acteristic lengths in the material parameters and/or in time are usually an exper- 
imentally robust feature. These scaling laws, and their characterizing exponents, 
are of interest to theoretical physics since they express a certain universality. At 
the same time, scaling laws (rather than more detailed features) are ameanable 
to heuristic and rigorous analysis and thus are a good test for the model and a 
challenge for mathematics. 

Scaling laws and their exponents reflect a scale invariance. In a multiscale 
model, these scale invariances are broken and only approximately valid in certain 
parameter and/or time regimes. The cross-over between two scaling laws reflects a 
change in the dominant physical mechanisms. In studying cross-overs, theoretical 
analysis may have an advantage over numerical simulation which has to explore 
many parameter decades and thus has to cope with widely separated length scales. 

Together with various collaborators, the author has analyzed scaling laws and 
their cross-overs in both static (variational) and dynamic models. The dynamic 
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models considered were of gradient-flow type and thus endowed with a variational 
interpretation: steepest descent in a multiscale energy landscape. The examples are 

• The branching of domains in uniaxial ferromagnets 1 (with R. Choksi and R. 
V. Kohn). Strongly uniaxial ferromagnets have only two favored magnetiza- 
tion directions ("up" and "down"). The width of the corresponding domains 
decreases towards a sample surface perpendicular to the favored axis. We rig- 
orously establish the scaling of the energy in the sample dimensions in support 
of this behavior. To leading order, the micromagnetic model behaves like a 
three-dimensional analogue of the Kohn-Muller |10| model for twin branching. 

• The period of cross-tie walls in ferromagnetic films [2] (with A. DeSimone, R. 
V. Kohn and S. Muller). Cross-tie walls are transition layers between domains 
in ferromagnetic films. They display a periodic structure in the tangential 
direction. The experimentally observed scaling of the period in the material 
parameters is not well- understood In this paper, we present a combination 
of heuristic and rigorous analysis which reproduces the experimental scaling 
and thus identifies the relevant mechanism. 

• The rate of capillarity-driven spreading of a thin droplet [S] (with L. Gia- 
comclli). Here, the starting point is the lubrication approximation. The scale 
invariant version of the model is ill-posed and has to be regularized near the 
contact line, c. g. through allowing finite slippage. In this paper, we rigorously 
derive a scaling law for the spreading of the droplet in an intermediate time 
regime. This scaling law depends only logarithmically on the length scale in- 
troduced by the regularization, in agreement with a conjecture of de Gennes 

• The rate of coarsening in spinodal decomposition (with R. V. Kohn). 
Spinodal decomposition is usually modelled by a Cahn-Hilliard equation. In 
the later stages, it is experimentally observed that the phase distribution 
coarsens in a statistically self-similar fashion. In this paper, we rigorously 
prove upper bounds for this coarsening process. The exponents are the ones 
heuristically expected and depend on whether the mobility is degenerate or 
non-degenerate: i 1 / 4 resp. t 1 ^. In PJ, we predict a cross-over for almost 
degenerate mobility due to a change in the coarsening mechanism. 

• The first-order correction to the Lifshitz-Slyozov- Wagner theory for Ostwald 
ripening [7] (with A. Honig and B. Niethammer). Ostwald ripening describes 
the late stage of spinodal decomposition in an off-critical mixture (volume 
fraction of one phase ^ < 1). The minority phase then consists of several 
particles immersed in a matrix of the majority phase. The particles are ap- 
proximately spherical and don't move — the Lifshitz-Slyozov — Wagner theory 
describes the evolution of the radii distribution. There is a major interest 
in identifying the next-order correction term in <fi. We rigorously show that 
there is a cross-over in the correction term from to (j) 1 / 2 depending on 
the system size. 

Our method to rigorously analyze these scaling laws in a multiscale model 
is based on relating integral quantities (energies, average length scales, dissipation 
rates...). It is different from the more local method of matched asymptotic expan- 
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sions. In particular, it differs from the latter by the absence of a specific Ansatz. In 
order to relate the integral quantities in our Ansatz-free approach, we need interpo- 
lation inequalities. These interpolation inequalities encode the competition of the 
dominant physical mechanisms in a scale-invariant fashion (e. g. the competition 
between driving energetics and limiting dissipation or between bulk and surface en- 
ergy). Hence tools from pure analysis are here employed in a more applied context. 
In order to illustrate this set of ideas, we present a simple application. 

2. An example from micromagnetics 

According to the well-accepted micromagnetic model, the experimentally ob- 
served ground-state of the magnetization m is the minimizer of a variational prob- 
lem. We are interested in transition layers ("walls") between domains in a film 
of thickness t in the {x\ 1 X2)-plane. We assume that the in-plane axis mi is fa- 
vored by the crystalline anisotropy so that domains of magnetization m = (0, 1,0) 
or m = (0,-1,0) form. In order to avoid "magnetic poles", the walls separating 
such domains are parallel to the ir^-axis. We are interested in their specific en- 
ergy per unit length in a^-direction. Hence the admissible magnetizations m are 
^-independent and connect the two end-states 

m = m(xi, X3) E S 2 for (x\. X3) G fi := (—00, 00) x (—5, |) 
and lim^^-too m 2 {x\, x 3 ) = ±1. 

The specific energy, which is to be minimized, is given by 

E{m) = d 2 [ \Vm\ 2 d 2 x + Q [ {m\ + m§) d 2 x + [ \Vu\ 2 d 2 x, (2.2) 
Jn Jn Jr 2 

where V refers to the variables x = (£1,2:3). Here the first term is the "exchange 
energy" , the second term comes from crystalline anisotropy and favors the m 2 -axis. 
The last term is the energy of the stray-field h s = — Vu determined by the static 
Maxwell equations 

V X h s = and V • (h s + m) = 0, 
which are conveniently expressed in variational form for the potential u 

[ m-V(d 2 x = / Vu-V(d 2 x for allC € C^°(J? 2 ). (2.3) 
Jn Jr 2 

We see that both "volume charges" ( V • m in Q) and "surface charges" (7713 on 
dfl) generate the field h s and thus are penalized. Since the energy density, i.e. 
|Vm| 2 , depends on m through (|2.3|l . the problem is non-local. The constraint of 
unit length, see l|2.1|l . makes the variational problem nonconvex. 

The model is already partially non-dimensionalized: The magnetization m and 
the field — Vu are dimensionless, but length is still dimensional. In particular, d has 
dimensions of length (the "exchange length" ) and Q is dimensionless (the "quality 
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factor"). Hence the model has two intrinsic length scales (material parameters), 
namely d and d/Q^, and one extrinsic length scale (sample geometry), namely t. 
Despite its simplicity, it is an example of a multiscale model and we expect different 
regimes depending on the two nondimensional parameters Q and 4. 

We will focus on the most interesting regime of "soft" materials (i. e. with low 
crystalline anisotropy) and thicknesses t close to the exchange length d 

Q < 1 and Q < (-) 2 < Q~ x . (2.4) 
d 

Numerical simulation suggest a cross-over within this range [SJ Chapter 3. 6, Fig. 
3.81]: 

• For thin films: "Neel walls" (see [5J Chapter 3.6 (C)]), whose geometry is 
asymptotically characterized by 

dm 

— — = and TO3 = =>■ m — (cos 9(xi), sin 9(xi), 0). (2.5) 
0x3 

• For thick films: "Asymmetric Bloch walls" (see $\ Chapter 3.6 (D)]), whose 
geometry is asymptotically characterized by 

— Vu = ==>■ V • m = in and 7713 = on dft 
=> ( mi ,m 3 ) = (-^,^) for a ^ with ^ = on 30. (2 ' 6) 

This cross-over in the wall geometry is reflected by a cross-over in the scaling of 
the specific wall energy E. Our proposition rigorously captures this cross-over in 
energy. 



Proposition 1 In the regime \2.4\j we have 



^m E{m) ~ W 2 ii /- W<^h \- (2J) 




By >, < we mean > resp. < up to a generic universal constant and ~ stands 
for both > and <. This scaling qualitatively agrees with the numerical study of the 

energy cross-over in the thickness 1 given in |51 Fig 3.79]. 

Upper bounds are proved by construction. Here we make the Ansatz l|2.5|) . 
resp. Q2.6[l. and let ourselves be inspired by the physics literature for the details of 
the construction. The matching lower bound in l|2.7|) states that one cannot beat the 
Ansatz — at least in terms of energy scaling — by relaxing the geometry assumptions 
(|2.5|l or 1)2. 6fl . Therefore Proposition ^ is a validation of the predicted cross-over in 
the geometry. We call this type of analysis Ansatz-free lower bounds. 

1 the x-axis corresponds to ^, the y-axis to J^, and Q = 0.00025 
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3. Proof 

The upper bound in Proposition ^ comes from the following two lemmas. We 
only sketch their proof since our main focus is on lower bounds. 

Lemma 1 For (^) 2 <C Q there exists an m of the form \2. b)) with 

E{m) ~ d 2 . (3.8) 

Lemma 2 For (^) 2 3> Q there exists an m of the form $2.5\) with 

E(m) ~ t 2 In" 1 J^. (3.9) 

For the lower bound we need to estimate the components mi and 7713 by E. 
In Lemma 01 we control 7713 by the stray-field and exchange energy. More precisely, 
the stray-field energy penalizes 7773 on <9f2 in a weak norm. We interpolate with the 
L 2 (£!)-control of Vm to obtain L 2 (fi)-control of 7713. In Lemma 0] we control the 
vertical average mi of mi by stray-field, exchange, and anisotropy energy. More 
precisely, the penalization of V • rn through the stray-field energy yields a penal- 
ization of in a weak norm. We interpolate with the L 2 (f2)-control of Vm 
(exchange) to obtain an estimate on the variation of mi . We then interpolate with 
the L 2 (fi)-control of mi (anisotropy) to obtain L°° (i?)-control of mi. 



Lemma 3 We have for any m satisfying \2.1\l 



E{m). (3.10) 



Lemma 4 In the regime (4) 2 3> Q we have for any m satisfying \2.1)) 



1 t 2 1 

sup mj( Xl ) < ( - In — + — ) E(m). i 3 . 1 I ) 



-oc,oo ) 



Proof of Lemma ^ The construction is due to Hubert [B]. We nondimen- 
sionalize length by t, i. e. t = 1. One can construct 2 a smooth tp; — > R with 

VtA| 2 < 1 inft, i> = on<90 and for \xx\ > 1, 

such that there exists a curve 7 C O with 

7 connects (0, — -|) to (0, |) and |VtA| 2 = 1 on 7. 

In line with the Ansatz l|2.6[) . we define m: fl — > S 2 via 

(-..-) = « = 1:1 f il W 



Ox*' 8*i " I + J I ri Stt 



2 Indeed, one possible recipe is to start from = h ~ \ x \ an< i to modify ij> outside of a 



neighborhood of the curve 7 = { ( J — &\ , %2 ) I #2 G [—51 5] } 
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Only exchange and anisotropy contribute to the energy: 

E{m) ~ d 2 + Q, 

which turns into (|3.8J) in the regime under consideration. 

Proof of Lemma [2] Making the Ansatz 12.5fl . the energy simplifies to 

f°° dm f°° f 

E(m) = d 2 t | \ 2 d Xl +Qt m\dx 1 + \Vu\ 2 d 2 x (3.12) 

J-oo dX! J_ QO J R 2 

^ <2 !t f r^trl^ + T I"" m 2 dx 1+ f \VU\ 2 d 2 x\ 1 
I t J^oo 1 - m{ dxi t J_ OQ J R 2 J 

where U is the harmonic extension 3 of mi from {x 3 = 0} onto J? 2 . Hence <|3~T2"1) 
holds for any extension U of mi. We now have to construct U such that its restric- 
tion mi satisfies m 2 (0) = 1 in order to allow for the sign change of m2- J R 2 | VZ7| 2 d 2 x 
just fails to control the L°°-norm of U and thus of mi — the counterexample involves 
a logarithm which we also use in this construction. The logarithm is cut off at the 
length scales \ <C <|: 



U(x) = In —g- lnymm{(— — ) 2 + (-^-) 2 , 1}. 

An elementary calculation shows l|3.9l) for mi(xi) = U(xi,0). A more detailed 
analysis of the reduced variational problem (|3.12l) is in [41 113). 
Proof of Lemma |3l We rewrite H2.3(l as 

d( 2 f du d( 2 / du d( 2 / 9mi 2 



m 3 -^z = / 7^7^ dx + ^J^ dx + ^Cd 2 x (3.13) 
/a 0x3 y fl 2 cte 3 dx 3 V fl 2 ctei ctei J a dx\ 

and choose the test function 

C(x!,X3) = m 3 (xi)r)(x a ) where x 3 = tx 3 

and ?7 £ Cq°(R) is chosen such that ^^(^3) = 1 for £3 6 (— |, |) in order to have 

9C/ \ 1 / s dr] 1 _ . . , t t 

— {xi,x 3 ) = -m 3 (xi) —{x 3 ) = ~m 3 (xi) tor x 3 e (--,-). 
0x3 r ax 3 r 2 2 

Hence the term on the 1. h. s. of (|3.13() turns into 

m 3 d 2 x = I m 2 dxi (3-14) 



dx 3 J_oo 
and the first term on the r. h. s. of l|3.13(l is estimated as follows 



du d( 2 
R2 dx 3 dx 3 



[\e I m\dxi ) . (:-!. K>) 



3 The inequality J fi2 |Vu| 2 ci 2 2 < t 2 / fi2 |V(7| 2 d 2 x can best be seen by expressing both integrals 
in terms of the Fourier transform mi(fci) of mi(ii). 
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The two remaining terms are also easily dominated 

du 



du d£ 
R 2 dx\ dxi 



d 2 x 



L 



dm 
R 2 dxi 



(d 2 x 



< 



< 



< 



< 



R i dx\ 
du 



,dm 3 s 



d^yd'xt c-^Ydx 



dxi 



(i?rd'x ( 



,dm 3 



. dm , , , 9 
(— ) 2 d 2 xt 
R 2 ax\ 



a dxi 



2 d 2 x 



< -E, (3.16) 

a 



mi dx\ 



t 



E I m 3 dx\ 



Collecting 13. 1411 - 1|3. 17(1 and using the Cauchy-Schwarz inequality gives 



11 t 

t ■ d + J 2 



rnldxi < ( y + -, + ) E < ( ^ + — ) E. 



1 t 
~t + d 2 



(3.17) 



(3.18) 



On the other hand, we use Poincare inequality in the ^-direction which we 
integrate over x% € (—00,00) 



dm 3 



t. 



(m 3 -m 3 ) 2 d 2 x < t 2 I {^Yd z x < {-) 2 E 



(3.19) 



Now (I3.18|) and l|3.19|l combine as desired into l|3.10|) . 

Proof of Lemma In the first step we establish for < p <C £ and 

< a - ii < t 

1 



6+p 



PJii 



mi dx\ — 



mi dx] 



11 1 
I - In - + — 

V p pt 



E. 



(3.20) 



In order to establish 1)3. 2U[) . we construct an appropriate test function £ for 12. 3|) . 
We first define C on the strip R x (— ~, |) as piecewise linear 

ii+p < Xi 

l^i-xi+p) ^1 < %i < fi+P 

C(xi,x 3 ) = { 1_ (1 < xi < £1 \. (3.21) 

\ {xi - lx + p) £1 - P < xi < ii 

X! < £1 - p 



£ is just defined such that 



m • VCd 2 .x = -- 



mi dxi H — 



mi c?a;i 



(3.22) 



p^i PJii- P 

For the r. h. s. of (|2.3(l we have to extend £ onto all of R 2 : We harmonically 
extend C on the upper and lower half-plane R x (|, +00) resp. 1? x (—00, — |). We 
claim 



/ \V(\ 2 d 2 x < In 



^ t 
- + -. 



(3.23) 
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This yields the following estimate of the r. h. s. of (|2.3I) 



/ Vit-V(d 2 j; 


< 


a 


J Ft 2 







Vu\ 2 d 2 x / \S7(\ 2 d 2 x 



< ( E (ml + t 
\ \ p p 

Obviously IlFE$ & j3~2ljl yields jOty . 

We now argue in favor of 13.23fl . On the strip J? x (— |, |) we have 



1^0 d x 2t/t>(-) 

p 



p 



(3.24) 



(3.25) 



The Dirichlet integral of the harmonic extension is estimated in terms of its bound- 
ary value as follows 



d x 



'flx(i, + oo) JO x 3 J-oo 

see [El Theoreme 9.4, Theoreme 10.2]. Since 



1 f°° t t 

{((xi+ x 3 ,-) - ((xi,-)) 2 dxidx 3 , 



t , 



(C(xi +x 3 ,-) ~((xi,-)) dx 



l3~2Tl 




this yields 



\V(\ 2 d 2 x < In-. 

7ix(i,+oo) ~ P 



Now (EH5|) and (|3~2rJj> combine into JOS- 

,2 

In the second step, we establish for t 3> %- and < £i — £i < 

1 , ^ 1 



For this, we observe that 
2 /-Si+p 



mi(xi) rfa?i - m\{i\) 



(3.26) 



(3.27) 



so that together with Q3.20|l we obtain 
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We now balance the first and last term by choosing p = ^- < I and so obtain (|3.27() . 
In the last step, we show <|3.11|) for -^-p 3> 1. For this we observe that 



f°° — 2 1 f 2 2 1 

/ m 1 dxi < — m l d x < —— E. 
J-oo ' in Qt 



Hence we obtain together with l|3.27[) for arbitrary £i E (-co, oo) 

mi{£,i) 2 < 7/ m 1 dxi + - (mi(6) - m 1 (x 1 )) 2 dx! 



,1 1 . It 1 . ^ 



Choosing € = <| ;§> ^- , we balance the two first terms and so obtain (|3.11(1 . 

Proof of Proposition^ It remains to establish the lower bound. For further 
reference we remark that by Poincare's inequality 

(mi -m l (0)) 2 d 2 x < t 2 [ |Vm 4 | 2 d 2 x < (^) 2 E. (3.28) 
~ Jn ~ d 

According to l|2.1|l . we have in particular \im Xl ^± 00 Wi2(xi) = ±1 and thus there 
exists an £i with m2(£i) = 0. W. 1. o. g. we assume £i = so that 7712(0) = 0. 
According to l|3.28|l we obtain 



L 



Furthermore, we have according to Lemma [3| 



m\d 2 x < i^fE. (3.29) 



/ m 2 d 2 x < (l + i-,) 2 ) E. (3.30) 



'(-I ,|)x(-| ,|) 

Since 1 - m\ = m 2 + m 2 ,, the estimates (I3.29jl & (I3.30jl imply 

/ (l-m 2 )d 2 x < + E. 

In view of l|3.28|) . this localizes to 

1 1 



l-m^Of <(- + -) E. (3.31) 



On the other hand, we have by Lemma0] 

1 h 

< 2 ^Qd? ' d? 
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provided (|) 2 3> Q. Combining l|3.31|l and l|3.32[) . we obtain 

/ 1 t 2 i i \ / l t 2 1 

1 < 72 ^TTIS + IS + 72 ) E ~ 7S^7T^ + 1S) E - (3-33) 



J 2 Qd 2 d 2 t 2 J \t 2 Qd 2 d 2 

Since we have by elementary calculus that 

T 2ln W 2 { > } h ln h{ > } ( ^ )2 ' 

l|3.33|) is equivalent to the lower bound in l|2.7|l . 

Acknowledgments. The author thanks A. DeSimone, Weinan E. , R. V. Kohn, 
and S. Muller for many stimulating discussions on micromagnetics. 
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